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Abstract. 

The dynamical behavior of a star network of spins, wherein each of N decoupled 
spins interact with a central spin through non uniform Hcisenbcrg XX interaction is 
exactly studied. The time-dependent Schrodinger equation of the spin system model is 
solved starting from an arbitrary initial state. The resulting solution is analyzed and 
briefly discussed. 
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1. Introduction 

In the last decades investigations on the properties of coupled-spin systems has gained 
an increasing interest in the quantum community|T]. Specially, the time evolution of a 
spin star system, that is a single spin coupled to a surrounding environment composed by 
a finite number of spins noninteracting [2]- [I] or interacting[5]-j6] among them, has been 
studied in detail. Central spin models provide an appropriate description of quantum 
information processes such as, for instance, quantum state transfer [7] and quantum 
cloning j8]. 

In this paper, we study in detail the dynamics of a XX central spin model that 
is composed by a localized spin 1/2, hereafter called central spin, coupled to an 
environment of N not interacting spins. The Hamiltonian describing our system is 

N N 

H = Luj2°i + «wi + E + (!) 

3=1 3=1 

The Pauli operators o~± refer to the central spin whereas the others, labelled by the 
index j, refer to the N environmental spins. The central spin Hamiltonian model 
([I]) can be successfully exploited to effectively describe many physical systems like 
quantum dots[9], two-dimensional electron gases jlO] and optical lattices |ll!j. The 
Hamiltonian model given by eq. (JTJ) is a realization of the so-called Gaudin model 
whose diagonalization has been derived in the framework of the Bethe ansatz (BA) 
[T2] . Such an approach provides however a rather formal solution whose implications 
for the dynamics of the spin system have not yet fully explored. In this paper we solve 
exactly the Schrodinger equation of motion of the total system starting from an arbitrary 
completely factorized initial condition. Our successfully treatment is strictly related to 
the circumstance that the component along the z axes of the total angular momentum 
operator S z — ^f- + | Z)j=i °~i = if" + Jz is a constant of motion. Our main result is the 
derivation of a closed formula for the time dependence of the probability amplitude of 
finding our spin system in any given state. 

2. Exact dynamics of XX-central spin model 

The goal of this Section is to derive the exact dynamics of our spin system starting 
from a completely factorized initial condition wherein the central spin, as well as p 
(p = 0, . . . , N — 1) of the N surrounding spins are in their respective up state, whereas 
the others ones are prepared in their down state. The case corresponding to p = N is 
trivial since the corresponding factorized state is an eigenstate of the Hamiltonian ([!])■ 

2.1. p=0 



The initial condition taken into consideration in this subsection is the following one: 

!#))) = IU)||...|>, (2) 
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where only the central spin is in the up state. It is easy to convince oneself [IJ that, 
thanks to the conservation of S z , the state ([2]) evolves into the state representable by the 



N 



following normalized superposition \ij)(t)) = a(t)\ 1a)\1 •••!)+ J2 
where 
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with A = — uoq. We underline that starting from such an initial condition, 
the time evolution is characterized by only one effective frequency, namely a e ff = 



J2j=i a 2 + A 2 . Thus the spin system fully restores its initial condition with a period 
T = 2ir/a e ff and behaves as if its dynamics were governed by an effective Hamiltonian 
model like that one given by eq. (JTJ), where Oij is substituted by a e //> independent of j. 



2.2. p=l,2. .., N-l 

Our aim is now to treat the more complicated dynamics of the XX central spin system 
starting from an arbitrary initial condition of the form 

\m) = \n)\i...n 1 ...n P ...i), (5) 

where p (p / 0) of the iV uncoupled spins, labelled by ki,...,k p , are in their up 
state | |), while the remaining N — p spins are in their down state | |). Since 
[S Z ,H] = 0, we claim that at any time instant t the system evolves into a normalized 

' N ' 



superposition of 



V 



= C N states wherein the central spin, as well as p among the 



N, are up and [ j = Cft 1 states wherein the central spin is down and p + 1 

V p + 1 / 

spins among the N are up. Thus, starting from the initial condition (jHJ), the vector 
state of the system evolves within a finite dimensional subspace whose dimension is 



+ c- 



■p+1 

N 



Cn+x- Starting from this property we proceed to write down effectively 



the evolved state of the system. To represent it, we exploit the set of p-tuples, 
that is the set of all the subsets of p elements from the first iV natural numbers, 



{(ii,i 2 , • • -,i P ), 1 < h < ... < i p 



< N}. It is well-known that the number of 
all p-tuples from iV numbers is exactly C%. Therefore, establishing a bijection between 
the set S p and the set of the states {| tx)|ltii • • • h P ■ ■ ■ !)}> as we ll as between the set 
S p+ i and the set {| 1a) II • • • Tji ■ ■ ■ ITjp+i)}) ^ is possible to represent the state of the 
system at the generic time instant t as follows: 



Sp+i 



32,—,3p+ 



oWIUUTix--- W,(6) 
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where the p-tuple (ii,i 2 , ■ ■ ■ ,i p ) G <S P identifies the probability amplitude ar ili i 2r „ ji \(t) 
of finding central spin A and exactly the spins ii,i 2 , . . . ,i p among the N around spins in 
their respective up state. Analogously the (p+l)-tuple (ji,j 2 , ■ ■ ■ ,jp+i) G <Sp+i provides 
the probability amplitude &(ji,j 2 ,..-,i P +i)(0 °f finding out the spin system in the particular 
state with the central spin down and exactly the environmental spins ji,j 2 , • • • , jp+i up. 
In order to get explicit equations for {a>(i 1) i 2 ,...,i p ){t)} and {b(j lt j 2 ,-,jv+i)(^)} we start from 
the time-dependent Schrodinger equation, introducing an appropriate mathematical 
notation useful to represent the transformations undergone by the states appearing 
in the expression ([6]) by the application of the Hamiltonian ([!])■ For this reason, we 
define two families of mappings {O r }^ =l and {S r }^ =1 . For any fixed r, the mapping O r 
transforms a p-tuple into a (p + l)-tuple accordingly to the rule 

O r : S p {r} -> Spl[\ O r (ii,i 2 , ■ ■ ■ ,i P ) = (ii,i 2 , • • • ,h) U {r}. (7) 

where «S~fr^(5pfr^) represents the set of all subsets of p elements, diverse by r (including 
r), from the the first iV naturals numbers. The mapping O r adds the natural number 
r to the p-tuple (ii, . . . , i p ), arranging them in increasing order. We point out that this 
correspondence is well defined if and only if r does not belong to the set {ii, . . . , i p }. 
The family of mappings {5 r }™ =1 on the contrary transforms, for any fixed r, a (p + 1)- 
tuple in a p-tuple in accordance with 

S r : Spl[ } S p ~ {r} , 5 r (ji,j 2 , ■ ■ ■ , jp+i) = (ji, 32, ■ ■ -Jp+l) - {r}- (8) 
It acts on the family of p + 1 elements, recovering a p-tuple from {ji, j 2 , . . . ,j p +i} 
by eliminating the element r. Obviously the above correspondence is well defined 
if and only if r belongs to the set {ji, j 2 , . . . , jp+i}- Inserting eq.® into the time- 
dependent Schrodinger equation, we obtain the following system of coupled equations 
for the probability amplitudes a(u,i 2) ...,ip)(^) an d ^0i,j 2 ,...jp+i)(0 

N 

ia { i lMr .. iip) {t) = Aa (l - 1>i2) ... )ip) (t) + a r b 0r{il ,..., lp ) (t) (9) 

r=l(r^{ii,i2,...,! p }) 
N 

ib Ui,h,-,j P +i)( t ) = - Ab (h,h,...,j P +i)(t) + a rasr(ji,...,j p+1 )(t). (10) 

r=l(r£{j 1 ,j 2 ,--;jp+i}) 

Solving the above system requires the diagonalization of the companion matrix of the 
system which is of order Cp^ 1 . On the other hand, we notice that the mean value of 
< > may be expressed as < Sf >= £ |a (ili ... iip )| 2 - 1/2 = 1/2 - £ |%i,...,i p+ i) \ 2 , 

so that to decouple the system of eqs. fl9l)- f|T0l is of physical and mathematical 
interest. Thus, to proceed further we follow a standard procedure by which we 
succeed in converting the above system into two decoupled systems for each unknown 
set {ci'(i 1 ,i2,...,ip){t)} and {b(j 1 j 2! ...j p+1 )(t)}. In this way we get two systems of coupled 
equations for the amplitudes a>(i lt i 2 ,...,i p )(t) and &(ji,j 2 ,-,i P +i)(*) respectively. We have 
indeed 

N 

'<(il,i2,-,i„)(X) + CtrCtsas^o^,...^))^) 
r,s=lr^s 
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'(ji>h,-,j P +i) 



At) 



a 2 + E 



p+1 \ N 

\ b (hJ2,...j P+ i)( t ) + Y. a rasb 0s {5 r {h,...,j p+1 )){t) 

i=l J r,s=lrj^s 



'12) 



Eq. lfTTI) (Eq. (|T2|) ) defines a linear system of C% (C^ +1 ) coupled second order differential 
equations in the variables o,(i 1 ,i 2 ,...,i p )('t) (^0'ij2,-,j P +i)(*))- The C^(C^r +1 ) amplitudes 
(ii,i2,...,ip)(*)' (%i,j 2 ..--,vfi) W) ma y °e or dered in accordance with lexicographical 
prescription, that is a^,...,^) (%;,.. follows a^,...,^) (6(i l5 ...,i p+1 ) ) if z'i = z'i,z 2 = 
i 2 , • • • , *m-i = C-i^m < i'm witn m = !> 2 > • • • >P ( m = 1,2, . . . ,p + 1). Therefore, 
eqs. (TTTT) . ( (fT2|) ) admits the matrix representation 

x(t) = -Xx{t) (13) 

where x(i) = a(t)(b(t)) is the lexicographically ordered vector of the probability 

amplitudes a(i 1( i 2 ,...,ip)(^), (^0'i,i2,...,i P +i)(*)) an d ^ = -4, (£>) is the corresponding 
companion matrix. In accord with eq. ( TTTT) . the matrix elements of A — 
(Anm')i< m ,m'<c p N are given by 



A 



A 2 



JV 

E a?, 



i/ card({m — m'}) 



0. 



z/ card({m — m'}) 
otherwise. 



0(<=> m 
1 



771 



(14) 



where m = (ii,i2, ■ ■ ■ ,i p ) and card({m — m'}) is the total number of the elements in 
the difference set {m — m'}. In a similar manner, denoting now the (p + l)-tuples by 
q = (ji, . . • , jp+i), the matrix B = (^ 9g ')i<g, 3 '<c p+1 is defmed by 



p+i 
A 2 + E « 2 



i=l 



Qq—q'Qq'—qi 



if card({q — q'}) = 0(<^ g = g') 
z/ card({q — q'}) = 1 



(15) 



0, otherwise. 

It is quite simple to notice that the matrices A and B are symmetric. For example, if 
p = 1 (p = N — 2) the matrix elements of A, (B) G Mjy [C] assume the following simple 
form 

A 



N 

A 2 + E« 2 •/ 



IJ 



J 



B 



a 



aiCtj 



if i^j 



N 

A 2 + E «, 2 

ajv+i-i^iv+i-j, 



AT+l-i) 
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J- 



*/ « ^ J- 



(16) 



Exploiting the well known solution of a matrix second order initial- value equation [T3] 
like eq. ([TBI yields: 



x(t) 



+oo / i \fc + 2fc 

\ ^ I -"-J r yfc 



x(0) + 



■+oo / i \fc + 2A;+l 
\ "* I l ) 1 yh 

to (2* + l)! 



x(0) 



Moreover, if A' is a non singular matrix, taking into account that X admits nonsingular 
square roots, the solution ( 1TT1) may be rewritten in the following closed form 

x(t) = cos (y/Xt) x(0) + sin (v^t) (v / ^) _1 x(0). (18) 

Practically, the possibility of exploiting the solution (fTTl) / (TT8]) depends of course on 
our ability of diagonalizing the matrices A and/or B. The dimensions C% and C^ 1 
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and the structure of the two matrices A and B respectively make anyway difficult an 
analytical treatment. From this point of view, solution (fTTl) might play a formal role 
only. However, it's always possible to recourse to the numerical diagonalization of the 
matrices A and/or B. We may wonder in this case about the efficiency of this numerical 
treatment in comparison with other numerical procedures[H]-[l5]. We observe however 
that such treatments get effective results when N is confined to values within 10, more or 
less, mainly due to exponentially increasing resources required from such computations. 
Overcoming these limitations on N becomes thus a mandatory target to improve the 
quality of the results achieved from a numerical approaches. The calculation scheme 
introduced by us is based on numerical diagonalization of the matrices A and B and 
presents the advantage that it enables the successful treatment of central spin models 
possessing a large number of bath spins, at least for small or large polarization, that is 
for small or large p. Indeed, in these cases the dimensions of the companion matrices 
A and B is such that the performance of numerical simulations is not obstructed by 
computational obstacles. We believe that our diagonalization procedure might be of 
help to investigate the behavior of our system as a function of other initial conditions. 

3. Conclusions 

Over the last few years, a whole variety of methods has been applied to study 
decoherence phenomena in the central-spin models. The appropriate version of the 
Bethe ansatz has allowed to integrate Hamiltonian model, but the results are rather 
formal and practically very difficult to handle. In this paper we have analyzed the exact 
dynamics of a central spin nonuniformly coupled through the Heisenberg interaction to a 
surrounding environment composed by iV spins. Considering arbitrary initial conditions 
we have determined a numerically manipulative general solution from which information 
about the full dynamics of our spin models may be extracted. 
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